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The transverse momentum (pt) dependence of J/ip production in heavy ion collisions is investi- 
gated in a transport model with both initial production and continuous regeneration of charmonia. 
The competition between the two production mechanisms results in a p t suppression in central 
collisions, the gluon multi-scattering in the initial stage leads to a high p t enhancement, and the 
regeneration populates J/ips at low pt region and induces a minimum in R,AA(pt)- These three 
phenomena are indeed observed in both 200 GeV Cu+Cu and Au+Au collisions at RHIC energy. 

PACS numbers: 25.75.-q, 12.38.Mh, 24.85,+p 

The J/ip suppression in relativistic heavy ion collisions is widely accepted as a signature of the formation of 
Quark-Gluon Plasma (QGP) at finite temperature and density PJ. Recently, the PHENIX and STAR collaboration 
presented the transverse momentum distribution of J/ip production in Au+Au 0, H[ and Cu+Cu fij collisions at 
RHIC energy ^/snn = 200 GeV. In particular the new data display three interesting features (see Figs|T]and[2]): 

a) The averaged transverse momentum square (p 2 ) is strongly suppressed in central collisions, which is very 
different from the SPS data || where (p 2 ) gets saturated for central Pb+Pb collisions at energy ^snn ~ 17.3 
GeV; 

b) The nuclear modification factor Raa is about 1/2 at low p t but becomes around unity at high p t ; 

c) There exists a minimum in Raa in low p t region. 

Can these phenomena tell us something new about the nature of J/ip production at RHIC energy? 

The p t broadening at SPS energy is generally attributed to the gluon multi-scattering in the initial state [|| 
and the leakage effect in the final state In a pA collision, a gluon of the proton scatters from target nucleons 
before it fuses with a gluon from the target to form a J /ip. The gluon rescattering in the initial state is treated 
as a random walk in transverse momentum and the observed (pf) is predicted to increase linearly with the mean 
length of the path of the incident gluon. In an AB collision, both gluons which fuse to a J/ip are affected by 
the rescattering. The leakage effect on J/ip production has already been considered 20 years ago [l|, [&, 9] . The 
anomalous suppression inside the QGP is not an instantaneous process, but takes a certain time. During this 
time the J/ips with high transverse momenta may leak out of the source of the anomalous suppression. As a 
consequence, low p t J/ips are absorbed preferentially but high p t J ftps can survive. 

While charm quark production at SPS energy is expected to be small, there are more than 10 cc pairs produced in 
a central Au+Au collision at RHIC energy and probably more than 200 pairs at LHC energy [l(| . The uncorrelated 
charm quarks in the QGP can be recombined to form J/iJjs. Obviously, the regeneration will enhance the J/ip 
yield and alter its momentum spectrum. The regeneration approach for J/ip production at RHIC has been widely 
discussed with different models, such as thermal production on the hadronizaton hypersurface according to statistic 
law pd| , the coalescence mechanism [l2| , and the kinetic model [l3T ] which considers continuous regeneration in a 
QGP. Recently, the J/r+transverse momentum distribution at RHIC energy was discussed [l4[ in the frame of a 
two-component model jl5| which includes both initial J/ip production through nucleon-nucleon (NN) interaction 
and regeneration. 

The medium created in high-energy nuclear collisions evolves dynamically. In order to extract information about 
the medium by analyzing the J/ip distribution, both the hot and dense medium and the J / ip production processes 
must be treated dynamically. Since the massive J/ips are unlikely fully thermalized with the medium, their phase 
space distribution is governed by a transport equation. In the transport approach [l6| which includes dissociation 
and regeneration processes in a QGP, the pt broadening due to initial gluon rescattering and the leakage effect in 
the final state are respectively taken into account through the initial condition and the free streaming term of the 
transport equation. To comprehensively treat the J/ip distribution in the transport approach, the J/ip transport 
equation is solved together with hydrodynamic equations which characterize the space-time evolution of the QGP. 
In this Letter, we investigate in the transport model the J/ip transverse momentum distributions in heavy ion 
collisions at RHIC and LHC energies, including the (pi) as a function of centrality and the Raa as a function of 
Pt- 

In pp collisions, the J/ips from the feed-down of ip' and \ c are respectively about 10% and 30% of the total final 
J/ips [17|. Since J /ip,ip' ,Xc) is heavy, the distribution function /^.(p t , x t , r|b) in central rapidity region and 
transverse phase space (pt, x t ) at time r and fixed impact parameter b is controlled by a classical Boltzmann-type 
transport equation (l6j 
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The second term on the left-hand side arises from free-streaming of \& with transverse velocity = p t /y/pf t i/t^, 
which leads to the leakage effect and is important for those high momentum charmonia. The suppression and 
regeneration in hot medium are reflected in the loss term ay (pt , x t , r | b) and gain term f3y (p t , x t , t| b) . Considering 
only the gluon dissociation process g + — > c + c, ay can be expressed as 

ay(p t , Xt ,r\h) = ^- J {2 ^f° E W$(s)f g ( Pg ,T, u)) 8 (T — T c ) / 9 (T* - T) , (2) 

where E g and -E* are the gluon and charmonium energies, f g = l/(e p s" M / T — 1) is the gluon thermal distribution, 
T(x t , r|b) and w(x t , r|b) are the local temperature and velocity of the hot medium, and is the transition prob- 
ability of the gluon dissociation as a function of s = (p + p g ) 2 calculated with perturbative Coulomb potential [li ]. 
For J/ip, the dissociation cross section reads 

with A a = (2 11 7r/27)(m^ej/^) -1 / 2 , where ui is the gluon energy in the rest frame of J/ip, m c is the charm quark 
mass, and cjm is the binding energy of J/ip- From what discussed in [ljj[, to take the relativistic effect into 
consideration and to avoid non-physics divergence in the regeneration cross section, we replace the J /ip binding 
energy by the gluon threshold energy. 

The step function in the numerator of the loss term a, controlled by the critical temperature T c of deconfmcment 
phase transition, means that the dissociation happens only in the QGP phase, namely a — for T < T c . Since 
the hadronic phase occurs later in the evolution of heavy ion collisions when the density of the system is lower 
compared to the early hot and dense period, we have neglected the hadronic dissociation. The step function in the 
denominator of a, characterized by the dissociation temperature T*, indicates that any \t cannot survive when 
the temperature of the QGP is higher than T*, namely a = oo for T > T*. From the lattice QCD simulation [20| . 
the J/ip spectral function remains almost the same when the temperature is between T c and T^, but the sharp 
peak of J /ip vanishes suddenly when the temperature reaches T^. The introduction of the dissociation temperature 
in a is also consistent with the idea of sequential charmonium suppression [l], [2l[ . 

The gain term (3 in the transport equation (JT|) , which is a function of c and c recombination transition probability 
Wj!^ , can be obtained from the lose term a using detailed balance [HI, [l(| . Instead of the gluon distribution f g 
in a, the charm quark distribution f c in (3 is assumed to be in the form of 



da 



dy 



T A (xt)T B (x t -b) 



/ c (Pc,x t ,r|b) = ^ ^w^w ^/ cp ( Pc |b), (4) 



V=0 T 



where da^ N /dy\ y= Q is the production cross section of charm quark pairs at central rapidity region in NN 
collisions, and Ta and Tg are the thickness functions of the two colliding nuclei defined as T(x t ) = 
lim^ 1 ^_ 00 . Z2 ^ 00 T(x t , z\, Z2) with T(xj, Z\,Z2) = fz 2 dzp(x t ,z) and p(r) being the Woods-Saxon nuclear density 
profile. From the PHENIX data, open charms carry an elliptic flow V2 ~ 0.1 [22j which means that charm quarks 
are likely thermalized with the medium, we take the charm quark momentum distribution f cp as the Fermi-Dirac 
function f cp (p c , T, u|b) - l/(e^/ T + 1). 

With the known suppression and regeneration terms a and /3, the transport equation can be solved analytically 
with the result [l6| 

h (Pt, X*. r|b) = h (Pt, x* - v*(r - to), r |b) e" A. «^a(p*,**-v(r-r'),r'|b) 

+ f dr'P (pt,x t - v*(r - t'), r'|b) e" ^ «*r"«(p*,xt-v(r-r"),r"|b) i (5) 

The first and second terms on the right-hand side indicate the contributions from the initial production and 
continuous regeneration, respectively. Both suffer anomalous suppression. The coordinate shift x t — > x t — v$At 
reflects the leakage effect during the time period At. 

Since the collision time for NN interactions at RHIC energy is about 0.1 fm/c and less than the starting time tq 
of the medium evolution which is about 0.5 fm/c, the nuclear absorption for the initially produced charmonia has 
ceased before the QGP evolution and is reflected in the transport process as the initial distribution fy (p t , x t , ro|b) 
of the solution ([5]). Considering a finite formation time of charmonia which is about 0.5 fm/c and larger than the 
collision time, the nuclear absorption can be safely neglected for the calculations in heavy ion collisions at high 
energy [23[. In this case, the initial distribution can be written as 

5 dc* f 1 / p 2 \~ 6 

fa (pt,x t ,r |b) = / dzAdz B pA(x-t, ZAjPB^-t ~ b, z B ) , ,, 1 + 77-fr — (6) 



with a normalized power-law momentum distribution extracted from the PHENIX data for pp collisions [24j , where 

7 NN I 



c?cr* N I 'dy\ v —o is the charmonium production cross section at central rapidity region in NN collisions, and (p1)n is 
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the averaged transverse momentum square after NN collisions [6| 

(j%} N (-xt,z A ,z B \h) = (pt)NN + a gN pQ 1 (T A (x.t,-oo,ZA) +T B (x t -b,z B ,oc)) (7) 

with p being normal nuclear density. For given values b and x t in the transverse plane, suppose a \f is produced 
at longitudinal coordinates za and zb in nuclei A and B, respectively. The two gluons which fuse to form the 'J 
carry transverse momentum from two sources: 1) Intrinsic p tl because they had been confined to nucleons. The 
intrinsic part is observable via NN — > \& process and leads to (p1)nn in Eq.((7]). 2) Collisional contribution to 
pt, because in a nuclear matter, the gluons traverse thickness Ta (xt, — oo, za) and Ts(xt — b, zb, oo) of nuclear 
matter in A and B, respectively, and acquire additional transverse momentum via gN collisions. This is the origin 
of the second term in Eq.([7]). The constant a g M is usually adjusted to the data for pA collisions. 

The local temperature T(xt,r|b) and fluid velocity M M (xt,r|b), which appear in the thermal gluon and charm 
quark distribution functions and control the suppression and regeneration region via the two step functions in a 
and /?, are determined by the (2 + 1) dimensional Bjorken's hydrodynamic equations for the medium evolution, 

d T E + V-M = -(E + P )/t , 

d T M x + V • (M x v) = -M x /t - d xP , 

d T My + V • (M y v) = -My It - d y p , 

d T R + V • (Rv) = -R/t (8) 

with the definitions E = (e + p)j 2 — p, M = (e + p)j 2 v and R = 771, where 7 is the Lorentz factor, and e,p and v 
are the local energy density, pressure and transverse velocity of QGP. 

To close the hydrodynamical equations we need to know the equation of state of the medium. We follow Ref.[2f| 
where the deconfined phase at high temperature is an ideal gas of massless u, d quarks, 150 MeV massed s quarks 
and gluons, and the hadron phase at low temperature is an ideal gas of all known hadrons and resonances with 
mass up to 2 GeV (26[. There is a first order phase transition between these two phases. In the mixed phase, the 
Maxwell construction is used. The mean field repulsion parameter and the bag parameter are chosen as K=450 
MeV fin 3 [25| and £? 1 / 4 =236 MeV to obtain the critical temperature T c = 165 MeV at vanishing baryon number 
density. The initial condition of the hydrodynamic equation at RHIC energy is the same as in Ref. |16| . 

Solving the local temperature and fluid velocity from the hydrodynamic equations and then substituting them 
into the transport solution (0) for J/ip,ip' and Xc we obtain the distribution function /j/^(pt, x t , 17 |b) for the final 
state Jhps at the freeze-out time 77. By employing the Cooper- Frye formula [2lJ with a longitudinal Hubble- like 
fluid [28j . any physical observable A for J/-0 can be estimated by integrating A ■ fj/^ over the freeze-out surface. 

We now calculate the nuclear modification factor Raa and averaged transverse momentum square (p 2 ) as 
functions of centrality decided by the number of participant nucleons N p for J/t/js produced in Cu+Cu and 
Au+Au collisions at RHIC energy. The following parameters are used in our calculations: the charm quark mass 
m c = 1.87 GeV (HI, the charmonium mass mj/^, =3.6 GeV and m,// = m Xa =3.7 GeV, the starting time of 
hydrodynamic evolution to=0.6 fm, and the charm quark and charmonium production cross section at \fs = 200 
GeV da% N /dy\ y =o = 120 pb and B u da% N / dy\ y=0 =26.4, 4.4, and 13.2 nb for * = J/ip, f and X c uM- 

The "J binding energy or the gluon threshold energy in hot and dense medium should be smaller than its value 
in vacuum. It is estimated to be less than 220 MeV in the QGP phase [3(J Hl|. We will take ej/^,=150 MeV. 
From our numerical results, a not very large deviation from this value does not lead to a sizeable change in Raa 
and {p 2 ) ■ The quasi-free cross section with a modified binding energy by the medium effect has been discussed in 
detail in Ref.hJ. Such a quasi-free process has a larger width than that of the gluon dissociation with binding 
energy in vacuum. 

The results of Raa(N p ) and (p 2 )(N p ) at central rapidity for Au+Au collisions are shown in FigQ]and compared 
with the RHIC data @, Q ■ The Raa is defined as the ratio of J / tps produced in a nuclear collision to that in 
a corresponding pp collision, normalized by the number of binary collisions. When the centrality increases, the 
contribution from the initial production drops down monotonously, due to the suppression in the QGP, while 
the number of regenerated J ftps goes up monotonously, because the number of charm quarks increases with the 
collision centrality. From the experimental data, there exists a flat region at 50< N p <170 with Raa — 0.6, and at 
N p > 170 the Raa continues to decrease with increasing centrality. This suppression structure can be reproduced 

in our transport approach by choosing the dissociation temperature Tf — T^ c = T c and Tj — 1.92 T c . After 

the plateau ends, the space-time region with temperature T > increases with centrality, and more J/ips arc 

fully eaten up by the extremely hot QGP. This is the reason of the further decrease of Raa at N p > 170. Since 
the QGP with T > Td exists only in the early stage of the fireball and lasts for a short time, the initially produced 
charmonia experience this stage, while most of the regenerated charmonia that are created later do not. This 
is why the Raa for regenerated J/ips is much less influenced by the dissociation temperature. The participant 
number fluctuations becomes important to Raa at extremely large N p , as discussed at SPS energy [321 ] . and may 
explain the deviation of our result from the data in very central collisions. 

We also calculated the Raa for Cu+Cu collisions and compared it with the RHIC data [33[. When the colliding 
energy is fixed, the temperature of the fireball is mainly controlled by the participant number, and therefore the 
result for Cu+Cu is almost the same as that for semi-central Au+Au collisions with N p < 110, where 110 is the 
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maximum participant number for a Cu+Cu collision. Since Cu is much lighter than Au, the fireball formed in 
Cu+Cu collisions cannot reach the temperature for full J/ip dissociation. 
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FIG. 1: The nuclear modification factor Raa and averaged transverse momentum square {pi) as functions of participant 
nucleon number N p in central rapidity region for Au+Au collisions at RHIC energy. Dot-dashed, dotted and solid lines 
represent respectively the calculations with only initial production, only regeneration and both contributions. The data are 
from PHENIX collaboration 



Most of the produced particles are with low momentum [3(, they dominate the centrality dependence of the 
yield, and the information carried by high p t charmonia is therefore screened in Raa{N p ). In order to understand 
the behavior of those high p t J/ips which are more sensitive to the production and suppression mechanisms, we 
consider the averaged transverse momentum square (pf) in the bottom panel of Fig[TJ In our calculation, the 
values of (p 2 t ) N N = 4.14 (GeV/c) 2 [H and a gN = 0.1 (GeV/c) 2 /fm Q are used for NN collisions, see eq.([7]). 
At SPS, the gluon rescattering parameter was taken as a g N = 0.077 (GeV/c) 2 /fm 0, [34j|. As one can see, the 
initial contribution to (pf) increases smoothly at N p < 170 and shows a saturation when N p > 170. This behavior 
is very similar to the case at SPS energy [7| where there is almost no regeneration and the initial production 
can be considered as the total result. The regenerated J ftps are from the thermalized charm quarks, therefore, 
their transverse momentum is rather small and the averaged value is almost centrality independent, in comparison 
with the initially produced charmonia, see the dotted line. Since the regeneration becomes more important in 
central collisions, the total (p 2 ) is strongly suppressed at large N p as a consequence of the competition between 
the initial production and regeneration. While almost all the models with and without regeneration mechanism 
can describe the J/ip yield after at least one parameter is adjusted, the pt suppression in central collisions seems 
to be a signature of charmonium regeneration in Au+Au collisions at RHIC energy. Again, the calculation for 
Cu+Cu collisions is very similar to the result for Au+Au in the region of N p < 110. 

To look into the details of J/ip production and suppression mechanisms, we now turn to the calculation of 
RAA(pt) as a function of pt which is discussed in different models [HI, The result for central Cu+Cu and 
Au+Au collisions at midrapidity is shown and compared with the RHIC data in Figj2] The monotonous increase 
of Raa with only initial production comes from three aspects. One is the pt dependence of the gluon dissociation 
cross section. The gluon energy relative to J/ip in the cross section (|3]) is scaled by J/ip binding energy which is 
much smaller than 1 GcV. Therefore, gluons with small energy are more likely to dissociate a J/ip, or in other 
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p t (GeV/c) 

FIG. 2: The nuclear modification factor Raa as a function of transverse momentum p t in central rapidity region for Cu+Cu 
and Au+Au collisions at RHIC energy. Dot-dashed, dotted and solid lines represent respectively the calculations with only 
initial production, only regeneration and both contributions. The data are from STAR for Cu+Cu and PHENIX [|| for 
Au+Au. 

words, J/ips with low momentum are easy to be eaten up by the hot medium. The second reason is the leakage 
effect with which the high momentum charmonia can escape the anomalous suppression region. The last one is 
the pt broadening happened in the initial state. Gluons get additional transverse momentum via collisions with 
participant nucleons before they fuse to J ftps. At SPS energy the J/ip Raa has been observed to exceed unity 
at high p t which can be well explained by the gluon multi-scattering Q. At RHIC energy, from our numerical 
calculation, the gluon rescattering is essential for the increasing of Raa in high p t region. 

In comparison with the initially produced J ftps which can carry high momentum from the hard process, the 
regenerated J/ips from thcrmalizcd charm quarks are distributed at low momentum region. In Figf^l these effects 
are clearly shown by dot-dashed and dotted lines for initial production and regeneration, respectively In the 
low pt region, the competition between the initial production which increases with p t and the regeneration which 
decreases with pt leads to a relatively flat structure for central Cu+Cu collisions. 

For the pt dependence of Raa for central Au+Au collisions, shown in the bottom panel of Figf^l while the 
trend of the curves are similar to that for Cu+Cu collisions, there are obvious features arisen from the stronger 
suppression and stronger regeneration in Au+Au collisions. In comparison with Cu+Cu collisions, the fireball 
formed in Au+Au collisions is much hotter and larger and lasts much longer, and the initially produced J ftps are 
strongly suppressed on their way out of the fireball. Since not all the regenerated J/rps are from the extremely 
hot medium, the effect of the stronger suppression in Au+Au collisions on the regeneration is not so important 
as for the initial production. Considering the fact that the regeneration is proportional to the square of the 
number of binary collisions, it is more important in central Au+Au collisions than in Cu+Cu collisions. At low 
Pt the regeneration contribution even exceeds the initial production. It is the competition between the stronger 
suppression of the initially produced J/ tps and the stronger regeneration at low p t that explains the overall reduced 
J/tp Raa and its minimum at p t ~ 2.5 GeV/c in Au+Au collisions. For central Cu+Cu collisions, there seems 
to have a minimum too, see the top panel of FigO However, due to the relatively small contribution from the 
regeneration, our calculations do not show such a clear minimum. 

As a prediction, we calculated the J/ip Raa and (pf) at LHC energy, the results are shown in Figf5J The input 
to the transport equation, namely the charmonia and charm quark production cross sections in NN collisions, is 
taken as da'^/dy = 2 fib and da^ N /dy = 0.7 mb, estimated by the CEM model and PYTHIA simulation [3(| . 
In comparison with the heavy ion collisions at RHIC, the formed fireball at LHC is hotter, bigger and longer lived, 
and almost all the initially produced charmonia are eaten up by the medium. On the other hand, there are much 
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FIG. 3: The nuclear modification factor Raa and averaged transverse momentum square (pt) as function of N p in cen- 
tral rapidity region for Pb+Pb collisions at LHC energy. Dot-dashed, dotted and solid lines represent respectively the 
calculations with only initial production, only regeneration and both contributions. 



more charm quarks generated at LHC energy, the regeneration controls the population in central and semi-central 
collisions. This regeneration dominance leads to an increasing Raa and a much stronger p t suppression at LHC 
for large enough N p . The rapid change of Raa and (pf) in the small N p region is due to the strong competition 
between the initial production and regeneration, and the degree of this competition decreases with increasing 
centrality and ends when N p is large enough. The increasing Raa and especially the small and saturated (jpf) for 
central collisions can be regarded as signatures of regeneration dominance at LHC, which are not expected in a 
model with initial production only. 

In summary, we investigated the transverse momentum dependence of J/ip production in heavy ion collisions 
at RHIC and LHC energies. By considering both initial production and regeneration and solving the coupled 
hydrodynamic equations for medium evolution and the transport equation for J/rp motion, we focused on the 
centrality dependence of (pl)(N p ) and pt dependence of RAA(pt) in Au+Au and Cu+Cu collisions. While the 
high p t behavior is characterized by the initial production, the regeneration and initial production become equally 
important at low p t . At RHIC, the competition between the two production mechanisms leads to the decrease of 
J/tp (pt)(N p ) in central collisions and the minimum of Raa{pi) at low p t . It is necessary to emphasize that the 
initial production alone cannot reproduce the p t suppression and the minimum of Raa found at RHIC. At LHC, 
almost all the initially produced charmonia are eaten up by the hotter, larger and longer lived fireball, the J/ip 
behavior is dominated by the regeneration when N p is large enough. As a result, the J/tp Raa increases with 
centrality and decrease with p t in central collisions. On the other hand, there is almost no regeneration at FAIR 
energy, and the J/tp production is governed by the initial production. In this case, both the pt suppression and 
the minimum of Raa will disappear in heavy ion collisions at FAIR energy. 
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